Epidemiological drivers of transmissibility and severity of SARS-CoV-2 in England (preprint)

As the SARS-CoV-2 pandemic progressed, distinct variants emerged and dominated in England. These variants, Wildtype, Alpha, Delta, and Omicron were characterized by variations in transmissibility and severity. We used a robust mathematical model and Bayesian inference framework to analyse epidemiological surveillance data from England. We quantified the impact of non-pharmaceutical interventions (NPIs), therapeutics, and vaccination on virus transmission and severity. Each successive variant had a higher intrinsic transmissibility. Omicron (BA.1) had the highest basic reproduction number at 8.1 (95% credible interval (CrI) 6.8-9.3). Varying levels of NPIs were crucial in controlling virus transmission until population immunity accumulated. Immune escape properties of Omicron decreased effective levels of protection in the population by a third. Furthermore, in contrast to previous studies, we found Alpha had the highest basic infection fatality ratio (2.8%, 95% CrI 2.3-3.2), followed by Delta (2.0%, 95% CrI 1.5-2.4), Wildtype (1.2%, 95% CrI 1.0-1.3), and Omicron (0.6%, 95% CrI 0.4-0.8). Our findings highlight the importance of continued surveillance. Long-term strategies for monitoring and maintaining effective immunity against SARS-CoV-2 are critical to inform the role of NPIs to effectively manage future variants with potentially higher intrinsic transmissibility and severe outcomes.

Genomic assessment of invasion dynamics of SARS-CoV-2 Omicron BA.1 (preprint)

SARS-CoV-2 variants of concern (VOCs) arise against the backdrop of increasingly heterogeneous human connectivity and population immunity. Through a large-scale phylodynamic analysis of 115,622 Omicron genomes, we identified >6,000 independent introductions of the antigenically distinct virus into England and reconstructed the dispersal history of resulting local transmission. Travel restrictions on southern Africa did not reduce BA.1 importation intensity as secondary hubs became major exporters. We explored potential drivers of BA.1 spread across England and discovered an early period during which viral lineage movements mainly occurred between larger cities, followed by a multi-focal spatial expansion shaped by shorter distance mobility patterns. We also found evidence that disease incidence impacted human commuting behaviours around major travel hubs. Our results offer a detailed characterisation of processes that drive the invasion of an emerging VOC across multiple spatial scales and provide unique insights on the interplay between disease spread and human mobility.

Lineage frequency time series reveal elevated levels of genetic drift in SARS-CoV-2 transmission in England (preprint)

Random genetic drift in the population-level dynamics of an infectious disease outbreak results from the randomness of inter-host transmission and the randomness of host recovery or death. The strength of genetic drift has been found to be high for SARS-CoV-2 due to superspreading, and this is expected to substantially impact the disease epidemiology and evolution. Noise that results from the measurement process, such as biases in data collection across time, geographical areas, etc., can potentially confound estimates of genetic drift as both processes contribute "noise" to the data. To address this challenge, we develop and validate a method to jointly infer genetic drift and measurement noise from time-series lineage frequency data. We apply this method to over 490,000 SARS-CoV-2 genomic sequences from England collected between March 2020 and December 2021 by the COVID-19 Genomics UK (COG-UK) consortium. We find that even after correcting for measurement noise, the strength of genetic drift is consistently, throughout time, higher than that expected from the observed number of COVID-19 positive individuals in England by 1 to 3 orders of magnitude. Corrections taking into account epidemiological dynamics (susceptible-infected-recovered or susceptible-exposed-infected-recovered models) do not explain the discrepancy. Moreover, the levels of genetic drift that we observe are higher than the estimated levels of superspreading found by modeling studies that incorporate data on actual contact statistics in England. We discuss how even in the absence of superspreading, high levels of genetic drift can be generated via community structure in the host contact network. Our results suggest that further investigations of heterogeneous host contact structure may be important for understanding the high levels of genetic drift observed for SARS-CoV-2 in England.

The origins and molecular evolution of SARS-CoV-2 lineage B.1.1.7 in the UK (preprint)

The first SARS-CoV-2 variant of concern (VOC) to be designated was lineage B.1.1.7, later labelled by the World Health Organisation (WHO) as Alpha. Originating in early Autumn but discovered in December 2020, it spread rapidly and caused large waves of infections worldwide. The Alpha variant is notable for being defined by a long ancestral phylogenetic branch with an increased evolutionary rate, along which only two sequences have been sampled. Alpha genomes comprise a well-supported monophyletic clade within which the evolutionary rate is more typical of SARS-CoV-2. The Alpha epidemic continued to grow despite the continued restrictions on social mixing across the UK, and the imposition of new restrictions, in particular the English national lockdown in November 2020. While from a case-number perspective these interventions succeeded in reducing the absolute number of cases of SARS-CoV-2 in the UK, the impact of these non-pharmaceutical interventions was predominantly to drive the decline of those SARS-CoV-2 lineages that preceded Alpha. We investigate the only two sampled sequences that fall on the branch ancestral to Alpha. We find that one is likely to be a true intermediate sequence, providing information about the order of mutational events that led to Alpha. We explore alternate hypotheses that can explain how Alpha acquired a large number of mutations yet remained largely unobserved in a region of high genomic surveillance: an under-sampled geographical location, a non-human animal population, or a chronically-infected individual. We conclude that the last hypothesis provides the best explanation of the observed behaviour and dynamics of the variant, although we find that the individual need not be immunocompromised, as persistently-infected immunocompetent hosts also display a higher within-host rate of evolution. Finally, we compare the ancestral branches and mutation profiles of other VOCs to each other, and identify that Delta appears to be an outlier both in terms of the genomic locations of its defining mutations, and its lack of rapid evolutionary rate on the ancestral branch. As new variants, such as Omicron, continue to evolve (potentially through similar mechanisms) it remains important to investigate the origins of other variants to identify ways to potentially disrupt their evolution and emergence.

Context-specific emergence and growth of the SARS-CoV-2 Delta variant (preprint)

The Delta variant of concern of SARS-CoV-2 has spread globally causing large outbreaks and resurgences of COVID-19 cases. The emergence of Delta in the UK occurred on the background of a heterogeneous landscape of immunity and relaxation of non-pharmaceutical interventions. Here we analyse 52,992 Delta genomes from England in combination with 93,649 global genomes to reconstruct the emergence of Delta, and quantify its introduction to and regional dissemination across England, in the context of changing travel and social restrictions. Through analysis of human movement, contact tracing, and virus genomic data, we find that the focus of geographic expansion of Delta shifted from India to a more global pattern in early May 2021. In England, Delta lineages were introduced >1,000 times and spread nationally as non-pharmaceutical interventions were relaxed. We find that hotel quarantine for travellers from India reduced onward transmission from importations; however the transmission chains that later dominated the Delta wave in England had been already seeded before restrictions were introduced. In England, increasing inter- regional travel drove Delta's nationwide dissemination, with some cities receiving >2,000 observable lineage introductions from other regions. Subsequently, increased levels of local population mixing, not the number of importations, was associated with faster relative growth of Delta. Among US states, we find that regions that previously experienced large waves also had faster Delta growth rates, and a model including interactions between immunity and human behaviour could accurately predict the rise of Delta there. Delta's invasion dynamics depended on fine scale spatial heterogeneity in immunity and contact patterns and our findings will inform optimal spatial interventions to reduce transmission of current and future VOCs such as Omicron.

Context-specific emergence and growth of the SARS-CoV-2 Delta variant (preprint)

The Delta variant of concern of SARS-CoV-2 has spread globally causing large outbreaks and resurgences of COVID-19 cases. The emergence of Delta in the UK occurred on the background of a heterogeneous landscape of immunity and relaxation of non-pharmaceutical interventions. Here we analyse 52,992 Delta genomes from England in combination with 93,649 global genomes to reconstruct the emergence of Delta, and quantify its introduction to and regional dissemination across England, in the context of changing travel and social restrictions. Through analysis of human movement, contact tracing, and virus genomic data, we find that the focus of geographic expansion of Delta shifted from India to a more global pattern in early May 2021. In England, Delta lineages were introduced >1,000 times and spread nationally as non-pharmaceutical interventions were relaxed. We find that hotel quarantine for travellers from India reduced onward transmission from importations; however the transmission chains that later dominated the Delta wave in England had been already seeded before restrictions were introduced. In England, increasing inter-regional travel drove Delta's nationwide dissemination, with some cities receiving >2,000 observable lineage introductions from other regions. Subsequently, increased levels of local population mixing, not the number of importations, was associated with faster relative growth of Delta. Among US states, we find that regions that previously experienced large waves also had faster Delta growth rates, and a model including interactions between immunity and human behaviour could accurately predict the rise of Delta there. Delta’s invasion dynamics depended on fine scale spatial heterogeneity in immunity and contact patterns and our findings will inform optimal spatial interventions to reduce transmission of current and future VOCs such as Omicron.

Genomic reconstruction of the SARS-CoV-2 epidemic across England from September 2020 to May 2021 (preprint)

Despite regional successes in controlling the SARS-CoV-2 pandemic, global cases have reached an all time high in April 2021 in part due to the evolution of more transmissible variants. Here we use the dense genomic surveillance generated by the COVID-19 Genomics UK Consortium to reconstruct the dynamics of 62 different lineages in each of 315 English local authorities between September 2020 and April 2021. This analysis reveals a series of sub-epidemics that peaked in the early autumn of 2020, followed by a singular jump in transmissibility of the B.1.1.7 lineage. B.1.1.7 grew when other lineages declined during the second national lockdown and regionally tiered restrictions between November and December 2020. A third more stringent national lockdown eventually suppressed B.1.1.7 and eliminated nearly all other lineages in early 2021. However, a series of variants (mostly containing the spike E484K mutation) defied these trends and persisted at moderately increasing proportions. Accounting for sustained introductions, however, indicates that their transmissibility is unlikely to exceed that of B.1.1.7. Finally, B.1.617.2 was repeatedly introduced to England and grew rapidly in April 2021, constituting approximately 40% of sampled COVID-19 genomes on May 15.

SARS-CoV-2 lineage dynamics in England from January to March 2021 inferred from representative community samples (preprint)

Genomic surveillance for SARS-CoV-2 lineages informs our understanding of possible future changes in transmissibility and vaccine efficacy. However, small changes in the frequency of one lineage over another are often difficult to interpret because surveillance samples are obtained from a variety of sources. Here, we describe lineage dynamics and phylogenetic relationships using sequences obtained from a random community sample who provided a throat and nose swab for rt-PCR during the first three months of 2021 as part of the REal-time Assessment of Community Transmission-1 (REACT-1) study. Overall, diversity decreased during the first quarter of 2021, with the B.1.1.7 lineage (first identified in Kent) predominant, driven by a 0.3 unit higher reproduction number over the prior wild type. During January, positive samples were more likely B.1.1.7 in younger and middle-aged adults (aged 18 to 54) than in other age groups. Although individuals infected with the B.1.1.7 lineage were no more likely to report one or more classic COVID-19 symptoms compared to those infected with wild type, they were more likely to be antibody positive 6 weeks after infection. Viral load was higher in B.1.1.7 infection as measured by cycle threshold (Ct) values, but did not account for the increased rate of testing positive for antibodies. The presence of infections with non-imported B.1.351 lineage (first identified in South Africa) during January, but not during February or March, suggests initial establishment in the community followed by fade-out. However, this occurred during a period of stringent social distancing and targeted public health interventions and does not immediately imply similar lineages could not become established in the future. Sequence data from representative community surveys such as REACT-1 can augment routine genomic surveillance.

An integrated analysis of contact tracing and genomics to assess the efficacy of travel restrictions on SARS-CoV-2 introduction and transmission in England from June to September, 2020 (preprint)

BackgroundMitigation of SARS-CoV-2 transmission from international travel is a priority. Travellers from countries with travel restrictions (closed travel-corridors) were required to quarantine for 14 days over Summer 2020 in England. We describe the genomic epidemiology of travel-related cases in England and evaluate the effectiveness of this travel policy. MethodsBetween 27/05/2020 and 13/09/2020, probable travel-related SARS-CoV-2 cases and their contacts were identified and combined with UK SARS-CoV-2 sequencing data. The epidemiology and demographics of cases was identified, and the number of contacts per case modelled using negative binomial regression to estimate the effect of travel restriction, and any variation by age, sex and calendar date. Unique travel-related SARS-CoV-2 genomes in the COG-UK dataset were identified to estimate the effect travel restrictions on cluster size generated from these. The Polecat Clustering Tool was used to identify a travel-related SARS-CoV-2 cluster of infection. Findings4,207 travel-related SARS-CoV-2 cases are identified. 51.2% (2155/4207) of cases reported travel to one of three countries; 21.0% (882) Greece, 16.3% (685) Croatia and 14.0% (589) Spain. Median number of contacts per case was 3 (IQR 1-5), and greatest for the 16-20 age-group (9.0, 95% C.I.=5.6-14.5), which saw the largest attenuation by travel restriction. Travel restriction was associated with a 40% (rate ratio=0.60, 95% C.I.=0.37-0.95) lower rate of contacts. 827/4207 (19.7%) of cases had high-quality SARS-CoV-2 genomes available. Fewer genomically-linked cases were observed for index cases related to countries with travel restrictions compared to cases from non-travel restriction countries (rate ratio=0.17, 95% C.I.=0.05-0.52). A large travel-related cluster dispersed across England is identified through genomics, confirmed with contact-tracing data. InterpretationThis study demonstrates the efficacy of travel restriction policy in reducing the onward transmission of imported cases. FundingWellcome Trust, Biotechnology and Biological Sciences Research Council, UK Research & Innovation, National Institute of Health Research, Wellcome Sanger Institute. RESEARCH IN CONTEXTO_ST_ABSEvidence before this studyC_ST_ABSWe searched PubMed, medRxiv, bioRxiv, Web of Science and Scopus for the terms (COVID-19 OR SARS-COV-2) AND (imported or importation) AND (sequenc* OR genom* or WGS). We filtered the 55 articles identified through this search and rejected any that did not undertake SARS-CoV-2 sequencing as part of an epidemiological investigation for importation into a different country. The remaining 20 papers were reviewed in greater detail to understand the patterns of importation and the methods used in each case. Added value of this studyThis is the first published study on importations of SARS-CoV-2 into England using genomics. Plessis et al., (2021) used a predictive model to infer the number of importations in to the UK from all SARS-CoV-2 genomes generated before 26th June 2020. The current study assesses the period 27/05/2020 to 13/09/2020 and presents findings of case-reported travel linked to genomic data. Two unpublished reports exist for Wales and Scotland, although only examine a comparatively small number of importations. Implications of all the available evidenceThis large-scale study has a number of findings that are pertinent to public health and of global significance, not available from prior evidence to our knowledge. The study demonstrates travel restrictions, through the implementation of travel-corridors, are effective in reducing the number of contacts per case based on observational data. Age has a significant effect on the number of contacts and this can be mitigated with travel restrictions. Analysis of divergent clusters indicates travel restrictions can reduce the number of onwards cases following a travel-associated case. Analysis of divergent clusters can allow for importations to be identified from genomics, as subsequently evidenced by cluster characteristics derived from contact tracing. The majority of importations of SARS-CoV-2 in England over Summer 2020 were from coastal European countries. The highest number of cases and onward contacts were from Greece, which was largely exempt from self-isolation requirements (bar some islands in September at the end of the study period). Systematic monitoring of imported SARS-CoV-2 cases would help refine implementation of travel restrictions. Finally, along with multiple studies, this study highlights the use of genomics to monitor and track importations of SARS-CoV-2 mutations of interest; this will be of particular use as the repertoire of clinically relevant SARS-CoV-2 variants expand over time and globally.

Environmental drivers of SARS-CoV-2 lineage B.1.1.7 transmission in England, October to December 2020 (preprint)

Previous work has shown that environment affects SARS-CoV-2 transmission, but it is unclear whether emerging strains show similar responses. Here we show that lineage B.1.1.7 spread with greater transmission in colder and more densely populated parts of England. We also find evidence of B.1.1.7's transmission advantage at warmer temperatures versus other strains, implying that spring conditions may facilitate B.1.1.7's invasion in Europe and across the Northern hemisphere, undermining the effectiveness of public health interventions.

Transmission of SARS-CoV-2 Lineage B.1.1.7 in England: Insights from linking epidemiological and genetic data (preprint)

The SARS-CoV-2 lineage B.1.1.7, now designated Variant of Concern 202012/01 (VOC) by Public Health England, originated in the UK in late Summer to early Autumn 2020. We examine epidemiological evidence for this VOC having a transmission advantage from several perspectives. First, whole genome sequence data collected from community-based diagnostic testing provides an indication of changing prevalence of different genetic variants through time. Phylodynamic modelling additionally indicates that genetic diversity of this lineage has changed in a manner consistent with exponential growth. Second, we find that changes in VOC frequency inferred from genetic data correspond closely to changes inferred by S-gene target failures (SGTF) in community-based diagnostic PCR testing. Third, we examine growth trends in SGTF and non-SGTF case numbers at local area level across England, and show that the VOC has higher transmissibility than non-VOC lineages, even if the VOC has a different latent period or generation time. Available SGTF data indicate a shift in the age composition of reported cases, with a larger share of under 20 year olds among reported VOC than non-VOC cases. Fourth, we assess the association of VOC frequency with independent estimates of the overall SARS-CoV-2 reproduction number through time. Finally, we fit a semi-mechanistic model directly to local VOC and non-VOC case incidence to estimate the reproduction numbers over time for each. There is a consensus among all analyses that the VOC has a substantial transmission advantage, with the estimated difference in reproduction numbers between VOC and non-VOC ranging between 0.4 and 0.7, and the ratio of reproduction numbers varying between 1.4 and 1.8. We note that these estimates of transmission advantage apply to a period where high levels of social distancing were in place in England; extrapolation to other transmission contexts therefore requires caution.

COVID-19 epidemic severity is associated with timing of non-pharmaceutical interventions (preprint)

Background: Unprecedented public health interventions including travel restrictions and national lockdowns have been implemented to stem the COVID-19 epidemic, but the effectiveness of non-pharmaceutical interventions is still debated. International comparisons are hampered by highly variable conditions under which epidemics spread and differences in the timing and scale of interventions. Cumulative COVID-19 morbidity and mortality are functions of both the rate of epidemic growth and the duration of uninhibited growth before interventions were implemented. Incomplete and sporadic testing during the early COVID-19 epidemic makes it difficult to identify how long SARS-CoV-2 was circulating in different places. SARS-CoV-2 genetic sequences can be analyzed to provide an estimate of both the time of epidemic origin and the rate of early epidemic growth in different settings. Methods: We carried out a phylogenetic analysis of more than 29,000 publicly available whole genome SARS-CoV-2 sequences from 57 locations to estimate the time that the epidemic originated in different places. These estimates were cross-referenced with dates of the most stringent interventions in each location as well as the number of cumulative COVID-19 deaths following maximum intervention. Phylodynamic methods were used to estimate the rate of early epidemic growth and proxy estimates of epidemic size. Findings: The time elapsed between epidemic origin and maximum intervention is strongly associated with different measures of epidemic severity and explains 46% of variance in numbers infected at time of maximum intervention. The reproduction number is independently associated with epidemic severity. In multivariable regression models, epidemic severity was not associated with census population size. The time elapsed between detection of initial COVID-19 cases to interventions was not associated with epidemic severity, indicating that many locations experienced long periods of cryptic transmission. Interpretation: Locations where strong non-pharmaceutical interventions were implemented earlier experienced much less severe COVID-19 morbidity and mortality during the period of study.

Genomic epidemiology of a densely sampled COVID19 outbreak in China (preprint)

Analysis of genetic sequence data from the pandemic SARS Coronavirus 2 can provide insights into epidemic origins, worldwide dispersal, and epidemiological history. With few exceptions, genomic epidemiological analysis has focused on geographically distributed data sets with few isolates in any given location. Here we report an analysis of 20 whole SARS-CoV 2 genomes from a single relatively small and geographically constrained outbreak in Weifang, People's Republic of China. Using Bayesian model-based phylodynamic methods, we estimate the reproduction number for the outbreak to be 1.99(95% CI:1.48-3.14). We further estimate the number of infections through time and compare these estimates to confirmed diagnoses by the Weifang Centers for Disease Control. We find that these estimates are consistent with reported cases and there is unlikely to be a large undiagnosed burden of infection over the period we studied.

Estimates of the severity of COVID-19 disease (preprint)

Background: A range of case fatality ratio (CFR) estimates for COVID 19 have been produced that differ substantially in magnitude. Methods: We used individual-case data from mainland China and cases detected outside mainland China to estimate the time between onset of symptoms and outcome (death or discharge from hospital). We next obtained age-stratified estimates of the CFR by relating the aggregate distribution of cases by dates of onset to the observed cumulative deaths in China, assuming a constant attack rate by age and adjusting for the demography of the population, and age and location-based under ascertainment. We additionally estimated the CFR from individual linelist data on 1,334 cases identified outside mainland China. We used data on the PCR prevalence in international residents repatriated from China at the end of January 2020 to obtain age-stratified estimates of the infection fatality ratio (IFR). Using data on age stratified severity in a subset of 3,665 cases from China, we estimated the proportion of infections that will likely require hospitalisation. Findings: We estimate the mean duration from onset-of-symptoms to death to be 17.8 days (95% credible interval, crI 16.9,19.2 days) and from onset-of-symptoms to hospital discharge to be 22.6 days (95% crI 21.1,24.4 days). We estimate a crude CFR of 3.67% (95% crI 3.56%,3.80%) in cases from mainland China. Adjusting for demography and under-ascertainment of milder cases in Wuhan relative to the rest of China, we obtain a best estimate of the CFR in China of 1.38% (95% crI 1.23%,1.53%) with substantially higher values in older ages. Our estimate of the CFR from international cases stratified by age (under 60 or 60 and above) are consistent with these estimates from China. We obtain an overall IFR estimate for China of 0.66% (0.39%,1.33%), again with an increasing profile with age. Interpretation: These early estimates give an indication of the fatality ratio across the spectrum of COVID-19 disease and demonstrate a strong age-gradient in risk.